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ABSTRACT 

Wc present phase resolved optical spectroscopy and Doppler tomography of 
V1341 Cygni, the optical counterpart to the neutron star low mass X-ray binary 
Cygnus X-2. We derive a radial velocity curve for the secondary star, finding a pro- 
jected radial velocity semi-amplitude of A'2 = 79±3 km s~^, leading to a mass function 
of 0.51 ± 0.06 Mq, ^30% lower than the previous estimate. We tentatively attribute 
the lower value of K2 (compared to that obtained by other authors) to variations in 
the X-ray irradiation of the secondary star at different epochs of observations. The 
limited phase coverage and/or longer timebase of previous observations may also con- 
tribute to the difference in K2- Our value for the mass function implies a primary mass 
of 1.5 ± 0.3 Mq, somewhat lower than previous dynamical estimates, but consistent 
with the value found by analysis of type-I X-ray bursts from this system. Our Doppler 
tomography of the broad He 11 A4686 line reveals that most of the emission from this 
line is produced on the irradiated face of the donor star, with little emission from the 
accretion disc. In contrast, the Doppler tomogram of the N ill A4640.64 Bowcn blend 
line shows bright emission from near the gas stream/accretion disc impact region, with 
fainter emission from the gas stream and secondary star. This is the first LMXB for 
which the Bowen blend is dominated by emission from the gas stream/accretion disc 
impact region, without comparable emission from the secondary star. This has impli- 
cations for the interpretation of Bowen blend Doppler tomograms of other LMXBs for 
which the ephemeris may not be accurately known. 
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1 INTRODUCTION 

Characterising the mass spectrum of neutron stars (NS) 
is critical to the understanding of the equation of 
state (EoS) which d escrib e s deg enerate nuclear matter 
iLattimer fc PrakashI |2004 120071 ). Low mass X-ray bi- 
naries (LMXBs) with NS primaries offer a promising 
means for increasing the number of such mass esti- 
mates, and better constrain i ng th e NS EoS. For example, 
Ivan den Heuvel fc Bitzarakil (|l995l ) have shown that some 
neutron stars could accrete up to 0.7 Mq over the course 
of the binary lifetimes. The detection of such massive neu- 
tron stars would significantly reduce the range of currently 
allowable NS equations of state. 

LMXBs themselves are sub-divided into transient and 
persistent systems. Persistent systems are those where the 
primary is continuously accreting at a significant fraction 
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of the Eddington limit, whereas in transient systems, oc- 
casional outbursts are separated by long periods of qui- 
esccnce, generally explai ned by the disc instability model 
(|Dubus. Hameurv fc Lasotai2001i 'l. 

In quiescent transient systems, the mass of the compact 
object can be constrained by dynamical studies if the sec- 
ondary is bright enough. Spectroscopic observations allow 
both the orbital period and the velocity of the secondary 
(projected onto the line of sight) to be determined, and 
these quantities are all that are required to calculate the 
mass function (/(M)), a lower limit to the mass of the pri- 
mary star. If the optical counterpart is sufficiently bright to 
allow higher resolution spectroscopy to be performed, the 
mass ratio can be estimated by measuring the rotational 
broadening of th e absorption lines originating in the sec- 
ondary star (e.g. iMarsh. Robinson fc Woodlll994l 'l. Finally, 
fitting of light curves can place constraints on the system's 
inclination, providing the final piece of information required 
to constrain the primary mass. 
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However, in LMXBs where the optical hght is domi- 
nated by the accretion disc, other techni ques are required to 
find th e system parameters. Specifically. ISteeghs fc Casarej 
found that in the case of Scorpius X-f , the Bowen 
blend emission (a blend of N ill and C ill emission lines near 
A4640) can be used to trace the motion of the secondary 
star. Doppler tomography of Bowen blend emission lines re- 
vealed a bright spot of emission, attributed to Bowen emis- 
sion from the irradiated surface of the secondary star. This 
technique has since been used to dete rmine the secondary 
velocity in a number of other systems jCasares et al.l [20041 : 
ICornelisse et al.ll2008l . and references therein), although we 
note that no independent verification of K2 has been ob- 
tained for any of these systems to date. 

C yg X-2 is a persis tent NS LMXB, discovered in 
1965 iBowver et al.l 1 19651 ). and is the second X-ray bi- 
nary f or which an optical c o unterpart (V1341 Cygni) was 
found (|Giacconi et al.lll967l ). [Orosz fc Kuulkeii (|l999l ) find 
that Cyg X-2 is located at a distance of ^7 kpc and 
the high X-ray flux implies that it is accreting at close 
to the Eddington limit. It is one of the few persistent 
systems where the se condary is visible, contr ibuting 70% 
of the optical flux llOrosz fc Kuulkerj Il999tl. Ahhough 
the s pectral type is A9 III ( Casares. Charles fc Kuulkerd 
19981). the mass of the secondary (M2) is ^0.6 Mq (see 
Podsiadlowski fc Rappapord |2000| . for more details). From 



the radial velocity (RV) curve base d on absorptio n line s pec- 
tra from the stellar photosphere, ICasares et al.1 l| 19981 ) es- 
tablished that Porb = 9.8444 ± 0.0003 d and K2 = 88.0 ± 
1.4 km s"'^, yielding a value for the mass function of the NS 
in this system of /(M) = 0.69 ± 0.03 Mq. They also mea- 
sured the systemic velocity (7) to be —209.6 ± 0.8 km s~^. 
Their estimate for the rotational broadening of the sec- 
ondary star of 34.2±2.5 km s~^ constrains the mass ratio, q, 
to be 0.34 ± 0.04. Setting the l ower limit of M2 in line with 
the pr edictions of iKing et al.l (|l997l ) allowed ICasares et al.l 
l| 19981 ) to constrain the primary mass (Mi) to be greater 
than 1.88 Mq, and the orbital inclination (i) to be less than 
61°, with 95% confidence. The existence of such a massive 
NS suggests that a stiffer EoS describes degenerate nuclear 
matter. 

Subsequently, by fitt ing U-, V- and B-band light curves, 
lOrosz fc Kuulkerg (|l999l ) found a value for the inclination 
of 62°. 5 ± 4°, which co mbined with the parameters found 
bv lCasares et al] (|l998t ) gave a value for the primary mass 
of 1.78 ± 0.23 Mq. This latter mass estimate is consistent 
with the canonical NS mass of 1.4 Mq at the 2a level 



iQrosz fc Kuulkerd (|l999l ) also found that X-ray heating of 
the secondary star contributed very little to the overall light 
curves, possibly because of the large orbital separation, or 
because a flared accretion disc is shielding the secondary 
from the X-rays. 

In contrast, analysis of X-ray bursts from Cyg X-2 al- 
lowed iTitarchuk fc Shaposhnikovl (|2002l ) to constrain the 
primary mass to be 1.44±0.06 Mq, lower than that obtained 
by previous dynamical studies, and equal to the canonical 
NS mass of 1.4 Mq, within the uncertainties. 

Here, we present results based on spectroscopic obser- 
vations of Cyg X-2 taken in 2006 September/October. Our 
primary motivation for this work was to perform Doppler to- 
mography for this long period system, to study the accretion 
flow, and especially to make the first comparison between 



the radial velocity of the secondary as inferred from the 
Bo wen blend technique and absorption line measurements 
fi.e. ICasares et a"l]|l998l ). Confirmation of this technique for 
a system where an independent measure of the radial veloc- 
ity is available gives further support for using Bowen blend 
to measure K2 in situations where the secondary star is not 
directly observable. In light of the NS mass estimates dis- 
cussed above, we also used our data to perform a new ra- 
dial velocity analysis of this system to confirm the suggested 
massive nature of the NS star primary. 



2 DATA 

Our data cons ist of optical spectra, a cquired with the FAST 
spectrograph jFabricant et al.lll998l ). at the Cassegrain fo- 
cus of the 1.5-m Tillinghast telescope at the Fred L. Whip- 
ple Observatory, Mt. Hopkins, Arizona. The FAST3 Uni- 
versity of Arizona STA520A (SN4377) CCD was used, with 
2720 X 161 15 \i.m pixels in the binned images. The spectra 
were taken using the 600 lines mm~^ grating, with a 2 arc- 
sec slit, giving a spectral coverage of ~2000 A. The spectra 
were obtained on 16 nights in 2006 September and Octo- 
ber, covering four orbits. One spectrum was obtained each 
night, along with HeNeAr arc lamp exposures, bias, flat and 
dark frames, and one spectrophotometric standard exposure 
(BD-l-28°4211). For 14 of the nights, only one arc spectrum 
was available for the wavelength calibration of each Cyg X-2 
spectrum. For the other two nights, the Cyg X-2 exposures 
were bracketed by arc spectra. 

The frames were firstly processed using the iRAlQ CCD- 
PROC routines to remove instrumental effects. The spectra 
were then optimally extracted using tasks in the iraf kp- 
NOSLIT package. Wavelength solutions were found by fitting 
a 3rd order cubic spline to ^^60 - 70 lines in each arc spec- 
trum, giving an RMS error of <0.04 A, and these solutions 
were then applied to the target frames. The resulting dis- 
persion was ~0.74 A pixeP^, with a resolution (measured 
from arc lines) o f ~2.3 A. For comparison, the observations 
of lCasares et al.l (|l998h had a resolution of ~0.5 - 0.8 A. Ta- 
ble[T]gives details of our observations. In order to verify the 
wavelength calibration, the sky lines were examined. Only 
five sky lines were isolated and strong enough to be useful 
i.e. with a FWHM similar to the spectral resolution. 

For nights where only a single arc spectrum was avail- 
able, we found that there were significant discrepancies be- 
tween the measured wavelengths of these sky lines and their 
rest wavelengths, as much as 0.5 A. For those nights where 
the target observations were bracketed by arcs, the inter- 
polated wavelength solution resulted in sky lines at their 
expected wavelengths. 

To corr ect for these offsets, we used the RVSAO pack- 
age in IRAF (|Kurtz fc Minkiri998l ) to cross-correlate the sky 
spectra against the sky spectrum from October 16 (which 
was calibrated using two arc spectra) and applied the mea- 
sured wavelength shifts to the target spectra. For the two 
nights where we had two arc spectra, we found that this 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. 16 normalized spectra from 2006 September and October, with a vertical offset of 0.4 units. The A9 III template spectrum is 
also shown. The wavelength regions masked from the cross-correlation are also marked. The spectra have been boxcar smoothed. 
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Table 1. Observations of Cyg X-2 from 2006 September and 
October. 



Date 


Exp. time 


Coverage 


Dispersion 


Orbital 


(ut) 




(A) 




(A) 


u b 

phase 


Sep 17 " 


1200 


3456 - 


5453 


0.743 


0.98 


Sep 18 " 


1200 


3455 - 


5453 


0.743 


0.08 


Sep 19 


1200 


3773 - 


5773 


0.744 


0.18 


Sep 20 


1200 


3780 - 


5779 


0.744 


0.28 


Sep 21 


869 


3779 - 


5780 


0.744 


0.38 


Sep 22 


1200 


3773 - 


5773 


0.744 


0.49 


Sep 25 


1200 


3770 - 


5770 


0.744 


0.78 


Sep 26 


1178 


3774 - 


5774 


0.744 


0.89 


Oct 03 


1200 


3775 - 


5775 


0.744 


0.60 


Oct 16 


1200 


3787 - 


5786 


0.744 


0.92 


Oct 17 


1200 


3786 - 


5785 


0.744 


0.02 


Oct 19 


1200 


3785 - 


5784 


0.744 


0.22 


Oct 20 


1200 


3788 - 


5788 


0.744 


0.33 


Oct 21 


1200 


3768 - 


5768 


0.744 


0.42 


Oct 22 


1200 


3780 - 


5780 


0.744 


0.52 


Oct 23 


1200 


3784 - 


5784 


0.744 


0.62 



" Data on September 17 and 18 were taken with a slightly 
different central wavelength to the other n i ghts. 

Based on the To fomid bv lCasares et all l[l99^ and Porb = 
9.84456 d. 



technique produced a wavelength scale very similar (maxi- 
mum difference of ~0.02 A) to that produced using the two 
arc spectra. 

After applying the sky line correction, the Cyg X-2 spec- 
tra were imported into the molly spectral analysis package, 
where the spectra were shifted to the heliocentric frame, 
outlying flux values were cleaned and the spectra were re- 
binned onto a common velocity scale of 50 km s~^ pixeP^. 
We found that the \a errors produced by iraf were sig- 
nificantly lower than the RMS of continuum regions, and 
that this was different for each of the spectra. Therefore, we 
increased the errors in each of the spectra to more closely 
match the RMS. The spectra were rectified by fitting the 
continuum with a spline function (after masking the major 
emission and absorption features), dividing by this fit and 
subtracting unity. Fig. [T] shows these spectra, each with a 
vertical offset of 0.4 units. 

We subsequently obtained spectra of 4 template stars 
in 2008 September , spanning the range A5 to Fl III (see 
ICasares et al] 119981 ). The A9 III template ( HR2489) is the 
same as that used by ICasares et al.l (|l998l ). These frames 
were processed in an identical manner to the Cyg X-2 
frames. However, arc spectra were taken at regular inter- 
vals, so the sky line correction procedure as employed for 
the Cyg X-2 spectra was not required. Inspection of the 
O III A5577.34 sky line showed that it was located at the 
correct wavelength in each calibrated frame. Details of the 
template observations are given in Table [5] The normalized 
spectrum for the A9 III template HR2489 is also shown in 

Fig.ni 

The spectra for each template were averaged and im- 
ported into MOLLY, where each spectrum was shifted to the 
heliocentric reference frame, and the continua removed as 
for the Cyg X-2 spectra. 



Table 2. Observations of template stars from 2008 September. 



Date 


Template 


Spectral 


No. 


Total exposure 


(ut) 


star 


type 


obs. 


time (s) 


Sep 24 


HD218260 


A5 III 


3 


180 


Sep 24 


HD220999 


A7 III 


4 


60 


Sep 24 


HD240431 


Fl III 


3 


270 


Sep 25 


HR2489 


A9 III 


9 


165 




4000 4500 
Wavelength (A) 



Figure 2. Flux calibrated spectrum from 2006 September 17 UT. 
The spectrum has been shifted to the rest frame of the system, 
and the major spectral features have been marked. 



3 RESULTS 

3.1 Flux calibrated spectrum 

For spectra where observing conditions were photometric, 
fiux calibration was performed against the standard star 
BD-f28°4211. Fig. [2] shows the fiux calibrated spectrum 
from 2006 September 17 from which we extracted a _B-band 
magnitude of 15.2 mag. This is consistent with the mean 
value over two decades reported bv lOrosz fc KuulkerEl(|l999l ) 
of 15.2 ± 0.2 mag (where the error refiects the RMS error of 
all the v alues). Also, visual c ompa rison of our spectrum with 
those of lvan Paradiis et al.l (|l990l ) and lO'Brien et al.l (|2004l ) 
shows that the spectrum is very similar at these epochs. 

The equivalent widths (EW) and full-width half- 
maxima (FWHM) of the absorption and emission lines in 
the average of all 16 spectra are shown in Table |3] 

The strongest feature in the Bowen blend in the av- 
erage spectrum is at a wavelength of A4637.4. Assumi ng a 
systemic velocity of —210 km s"'^ (|Casares et al.|[l998l ). we 
identify this emission as from the N ill A4640.64 line. This 
is freq uently observed to b e the strongest line in the Bowen 
blend (|Casares et al.|[2"004h . There does not appear to be any 
significant emission from the other lines normally present 
in the Bowen blend. However, because of the resolution of 
our spectra, emission from these lines may be present but 
blended with that from the N ill A4640.64 line. 
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Table 3. EW and FWHM of the emission and absorption lines 
in the averaged Cyg X-2 spectrum. 



Line 


EW 


FWHM 




(A) 


(A) 


He 11 A4686 


-2.9 ±0.1 


6.4 ±0.2 


Bowen blend 


-0.8 ±0.1 




H/3 


5.2 ±0.1 


33 ± 2 


H7 


4.4 ±0.1 


18 ± 2 


H5 


4.5 ±0.1 


18 ± 2 


HC 


3.6 ±0.1 


18 ± 2 


Ht? 


6.1 ± 0.1 


19 ± 2 



3.2 Radial velocity study 

As the primary motivation for this work was to produce 
Doppler tomograms, we initially did not observe any ra- 
dial velocity templates. However, by using the individual 
spectra as templates, we were able to derive a radial veloc- 
ity curve, which suggested that K2 was significa ntly lower 
than the value obtained by ICasares et al. I 1 19981 1. In order 
to confirm this lower value of K2, we later obtained spectra 
of 4 template stars, including the A9 III template used by 
ICasares et al.l (|l998l ). and we cross-correlated the Cyg X-2 
spectra against these templates to obtain radial velocities. 

For the cross-correlation, the broad Balmer absorption 
features, Ca 11 H and K lines, and emission lines were masked 
(see Fig. The trailed spectrogram in Fig. [3] shows a 
subset of the data used in the cross-correlation (AA4900 - 
5000). The S- waves due to the narrow absorption lines from 
the secondary are clearly visible. The cross-correlations also 
took into account the uncertainties in the wavelength so- 
lution of each of the spectra, as well as the slight differ- 
ences in the cross-correlation value when different interpola- 
tions/binnings were used. This was accomplished by running 
10000 cross-correlations for each spectrum-template pair. 
For each cross-correlation, the wavelength scale of the spec- 
trum and template were shifted independently by choosing 
a shift value from a Gaussian distribution with a mean of 
and standard deviation equal to the estimated uncertainty 
in the wavelength solution due to the arc line fitting and the 
sky line correction. 

We fit the 4 RV curves (derived from cross-correlation 
against the 4 templat es) using both circular a nd elliptical or- 
bit models. Using the lLucv fc SweenevI ()l97l[ ) test, we found 
that an elliptical orbit was not a significant improvement 
over a circular orbit for any of these 4 RV curves: the ec- 
centricity (e) in each case was ~0.06±0.04. With a circular 
orbit and 3 free para meters (if 2, 7 a nd Pp rbi with To set 
to the value found by ICasares et all (|l998l l). a sine fit to 
each of these 4 RV curves gave almost identical results (to 
bet ter than la). Using the same A9 III template as used 
by ICasares et all (|l998h . we found /C2 = 79 ± 3 km s~\ 
Porb = 9.84456 ± 0.00012 d and 7 = -212 ± 2 km s-\ 
where the radial v elocity of the template of 7.7±0.1 km s~^ 
(|Gontcharovll2006l ) has been included. This is in contrast to 
the values found by ICasares et al. I (|l998t ). For comparison, 
ICasares et al] (|l99^ found values of K2 = 88.0± 1.4 km s"^ 
Porb = 9.8444 ±0.0003 d and 7 = -209.6 ±0.8 km s'V The 
uncertainties in our measurements, here and in the remain- 
der of the paper, are Icr and were determined using Monte 
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Figure 3. Trailed spectrogram for the data between AA4900 and 
5000, showing some of the narrow absorption lines from the sec- 
ondary. 



Carlo simulations assuming Gaussian statistics. The values 
we obtained for Porb and 7 are consistent at approximately 
the Ic r level with the values determined by ICasares et al.l 
(| 19981 '). However, the K2 values differ at the 99.7% level. The 
RV curve obtained from cross-correlation against the A9 III 
template star and best-fitting circular orbit are shown in 
Fig. 131 The reduced {xi) for this fit is 1.2 for 13 degrees 
of freedom. For c omparison, the RV c urve generated using 
the parameters of ICasares et al.l l|l998l ) is also shown. 

3.3 System parameters 

Using the mass function equation and the values for K2 and 
Porb gives a value for the mass function of f{M) — 0.51 ± 
0.06 M|^. Because o f the higher resolution spectroscopy of 
ICasares et al.l (|l998l ). we use their value for the rotational 
broadening of the secondary star, vs'mi = 34.2±2.5 km s~^. 
Combined with our lower value of K2, we find a mass ratio of 
q = 0.42 ±0.06. Along with our value for f(M), we calculate 
Misin^i = 1.03 ±0.15 Mq. 

X-ray eclipses have not been observed from this system, 
and this allows us to set a firm upper limi t to th e inclination. 
Re-arranging the equation of iPaczvriskil (|l97ll ). we get: 

0.462 (^)^ 

where a is the binary separation and r2 is the radius of the 
secondary star. 

In order for the secondary not to eclipse the primary, 
the upper limit on the inclination is 90° — tan~^ which is 
72°. 7 ± 0°.6. This angle is relatively insensitive to changes 
in q. Using this limit, we place a lower limit on the primary 
mass of Ml > 1.18 ± 0. 17 Mp . 

lOrosz fc Kuuito^ 119991) find a lower limit to the or- 
bital inclination by fitting the ellipsoidal modulation in their 
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Orbital phase 

Figure 4. Radial velocity curve derived using the A9 III star 
HR2489 as a template. The best-fitting circular orbit is also plot- 
ted, and the 7 va lue of —212 km s~^ is marked with a dashed 
line. The curve of ICasares et al is plotted with a dotted 

line. Two phases are shown for clarity. 



l^-band light curve, assuming no contribution from the disc 
or from an irradiated secondary - both of these would di- 
lute the observed modulation, and so the inclination found 
when these effects are not inc luded is a lower limit. Ta king 
the published y-band data of lOrosz fc Kuulkerd (Il999l). we 
used the elc (Eclipsing Light Curve: Qrosz fc Hauschildtl 
cod e to fit the modulat i on, u sing the parameters de- 
scribed in lOrosz fc Kuulke"r3 (|l999|) , and reproduced their 
value for the inclination lower limilQ of ~49° (48°.7± 2°.l). 
We then adjusted the ELC parameters for Porb, .f{M) and 
q based on our spectroscopic results, and found the best fit 
was for i = 49°. 8 ± 2°.0. This sets an upper primary mass 
limit of 2.3 ±0.4 Mq. 

We then re-fit both the B- and \/-band light curves of 
lOrosz fc Kuulke"r3 |l99^ with ELC, with X-ray irradiation 
and the contribution from the accretion disc included. We 
adjusted the parameters for Porb, f{M) and q as before, and 
by optimizing the inclination and disc temperature param- 
eterqj found a best-fitting incli nation of z = 63 ± 3° , alm ost 
identical to the value found bv lOrosz fc Kuulkerd (|l999l l of 
i = 62.5 ± 4°. The uncertain t y was computed in a similar 
manner to lOrosz fc KuulkersI '(|l999l ). Combining this incli- 
nation estimate with our value of Mi sin"^ i — 1.03±0.15 Mq, 
we obtained a primary mass of M\ = 1.5 ± 0.3 Mq, and a 



secondary mass of M2 = 0.63 ± 0.16 M, 



0- 



2 Note that lOrosz fc KuulkersI l ll999l) found lower limits of 42° 
and 49° for the B and V-band light curves, respectively, and so 
adopted 49° as a lower limit. 

^ The equation in ELC which describes the dis c temperature as a 
function of radius is different to that used bv lOrosz fc KuulkersI 
1I1999I '). so we therefore firstly adjusted the temperature pa- 
rameters to give the sam e disc temperature profile as used by 
lOrosz fc KuulkersI |l99^ before optimizing these parameters. 



3.4 Doppler tomography 

Doppler tomography is an inversion technique where a ve- 
locity space image of the accretion disc, gas stream and sec- 
ondary star is g enerated using phase resolved spectroscopy 
(|Marsh fc Hornc 1988,; Marsh 20oT| ); the observed spectra 
are projections of the system at different orbital phases. A 
Doppler tomogram is a two dimensional image in velocity 
space (Kc,Vy), with intensity representing strength of emis- 
sion at a particular velocity. The Doppler tomograms pre- 
sented here were constructed using the maximum entropy 
method (MEM), as implemented in the doppler package. 
Synthetic line profiles, created from an initial tomogram of 
constant pixel values, were fitted to the observed line pro- 
files, by varying the pixel values. This process was contin- 
ued until the of the fit approached a target value. As 
a large number of tomograms can satisfy a particular tar- 
get the tomogram with the maximum entropy (i.e. the 
smoothest image) was chosen. The entropy was defined rel- 
ative to a Gaussian blurred default image. We obtained sim- 
ilar tomograms using the filtered back-projection technique 
implemented in molly, where the tomogram is created by 
smearing the filtered line profiles across the image in direc- 
tions determined by the orbital phases of the profiles. 



3.4-1 Hen X4686 emission line 

Fig. [S] shows the Doppler tomogram for the He 11 A4686 
emission line. This emission line is often seen in X-ray bina- 
ries, and is formed by reprocessing of soft X-rays from the 
central source. The Roche lobe of the secondary is plotted, 
for K2 = 79 km s~^ and q = 0.42 (see Section [3. 2|l . as is the 
centre of the mass of the system, marked by a *. The lower 
curve is the gas stream velocity from the inner Lagrangian 
point, and the upper curve is the velocity of the accretion 
disc along the gas stream. The main feature to note is that 
most of the emission comes from near the secondary star, 
although there is also a much broader region of faint emis- 
sion, due to the accretion disc, with a distribution typical of 
that seen in accretion discs in other LMXBs. 

Fig. [7] also shows the Doppler tomogram for the 
He 11 A4686 emission line, focussing on the emission from 
near the secondary star - the colour bar has been scaled 
to show only the peak emission. This tomogram shows that 
the peak of the emission arises on the inner face of the sec- 
ondary star. The centroid of this emission is at a velocity of 
Ky = 64 ± 15 km s"\ 



3.4-2 Nni \464O.64 emission line 

Fig. [8] shows the trailed spectrogram for the region near 
A4640. Because the systemic velocity has not been shifted 
out, the N III A4640.64 line appears near A4638. Although 
the S/N is low, it does appear that the S-wave for this line 
leads some of the other S- waves present, (which are due to 
the motion of the secondary star). 

Fig. [9] shows the corresponding Doppler tomogram, 
again with the Roche lobe and gas stream overplotted. 
The strongest emission is from the location of the gas 
stream/accretion disc impact region, with broad, fainter 
emission from the gas stream and near the secondary. This 
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Figure 5. Trailed spectrogram of the data used to construct the 
He II Doppler tomograms. 
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Figure 7. Doppler tomogram of the He II A4686 emission line, 
focusscd on the region of the secondary star. The Roche lobe of 
the secondary is overplotted for K2 = 79 km s~^ and q = 0.42. 
The position of the secondary is marked with an error bar. 



is consistent with the situation we see in the trailed spectro- 
gram (Fig. [8]). The overplotted S-wave has a phase offset of 
45°, similar to what we see in the Doppler tomogram. 
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Figure 6. Doppler tomogram of the He II A4686 emission line. 
The gas stream and Roche lobe of the secondary are overplotted 
for K2 = 79 km s'^ and q = 0.42. 



4 DISCUSSION 

4.1 Projected velocity of the secondary star 

Our value for K2 differs from the mo st recent observations 
of this system bv lCasares et aL I (|l998t ): Table H summarises 
the values f or 7 and K2 found by various authors. Because 
the work of ICasares et al.l l|l998h and our work show that 
the orbit is most likely circular, we also re- fi t the data 
published by Cowley. Crampton fc Hutchingj (| 19791 ) and 
ICrampton fc CowlevI (|l980l ) assum ing a circular orbit. 

In our work, and in the work of lCasares et al. I (|l998t ) the 
non-metal lines and the Ca 11 lines were masked out during 
cross-correlation (see Fig.[l|. Therefore all values of K2 and 
7 in Table 3] are essentially based on the velocities of the 
metal absorption lines. 

We now examine po ssible reasons for our lower value of 
K2 compared to that of lCasares et al.l (|l998l ). Firstly, it is 
possible that there are systematic errors in our wavelength 
solution, given that for most of our spectra only a single 
arc lamp exposure was available, and a correction based on 
the sky spectra was required. However, our RV curves ap- 
pear to be quite sinusoidal, which would be unlikely if each 
spectrum had random errors in the wavelength calibration - 
indeed, the RV curve from the uncorrected spectra deviates 
considerably from a sinusoid. Also as we mention in Section 
[21 for the cases where we had two arc spectra, the maximum 
difference between the wavelength solutions using the two 
arc spectra, and using the sky line correction technique is 
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Table 4. K2 and 7 found by different authors. 



Source 


Start of 


End of 


Time base 


Method used 


K2 


7 




observations 


observations 






(km s^-"-) 


(km s-i) 


a 


1975 June 03 


1978 August 09 


3.2 yr 


Metal lines 


87 ±3 


-222 ± 2 


b 


1975 June 03 


1978 August 09 


3.2 yr 


Metal lines 


89 ± 5 


-222 ±4 


c 


1979 July 02 


1979 September 17 


77 d 


Metal lines 


80 ± 6 


-203 ± 3 


d 


1979 July 02 


1979 September 17 


77 d 


Metal lines 


75 ±4 


-202 ± 2 


e 


1993 December 16 


1997 August 07 


3.6 yr 


Cross correlation 


88.0 ± 1.6 


-209.6 ±0.8 


f 


2006 September 17 


2006 October 23 


35 d 


Cross correlation 


79 ± 3 


-212 ± 2 



(a) KPNO data of ICowlev et al.l l|l979l '). 

fb^ Based on our re- fittin g of t he KPNO data of lCowlev et al.l lll979t) assuming a circular orbit. 
(cl lCrampton fc Cowled lll980l') . 

fd) Based on our re-fitt ing of the data of lCrampton fc CowlevI (|l980h assuming a circular orbit. 

(e) ICasares et al] lll998l) . 

(f) This paper. 
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Figure 8. Trailed spectrogram of the data near the N III A4640.64 
emission line. The systemic velocity has not been removed, so the 
N III A4640.64 line is near A4638. The solid line marks the S-wavc 
from N III A4640.64 emission, while the dotted line marks one of 
the nearby secondary absorption lines. 
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Figure 9. Doppler tomogram of the N III A4640.64 emission line. 
The gas stream and Roche lobe of the secondary are overplotted 
for K2 = 79 km s'^ and q = 0.42. 



~0.02 A, or ~1.3 km s"^ at A4686. Our RV curves include 
the uncertainties in wavelength scales, and also the slight 
differences in the cross-correlation values that arise when 
different interpolations are used. Fina lly, we find a value fo r 
7 almost identical to that found by ICasares et al.l (|1998D . 
which would again be improbable if there were any system- 
atic errors in our wavelength scales. 

Hence, we believe our measurement of K2 to be accu- 
rate, and so we require an alternative explanat ion for the dif- 
ference between this measurement and that o f[c asares et al.l 
(|l998h . Because our He 11 Doppler tomography confirms 
that the secondary is irradiated by X-rays from the cen- 
tral source, the absorption lines from the irradiated face of 
the secondary star may be weaker than those from the non- 



irradiated face. The irradiation of the secondary may be 
such that the absorption lines are completely quenched on 
the irradiated hemisphere of the secondary star. In this case, 
an RV curve based on the velocity of absorption lines from 
the secondary star will result in a K2 value which reflects 
the velocity of the non-irradiated hemisphere. This is larger 
than the velocity of the centre of mass of the secondary star. 
If the irradiation of the secondary was variable, such that 
at some phases the absorption lines on the irradiated face of 
the secondary were not fully quenched, then one would ex- 
pect the K2 value derived from the resulting RV curve to lie 
somewhere between the true value a nd the maximum value 
obtained for complete quenching. See Wade fc Hornd 
for further discussion of this effect. [Steeghs fc Jonkerl 



198i 



200 



argue that differences in K2 found by different authors for 
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the LMXB 2S 0921—630 may be due to variable irradiation 
effects, and we s uggest that the same m ay well apply in the 
case of Cyg X-2. IClarkson et aD (|2003l l show that the long 
term X-ray variability in Cyg X-2 can be explained by a 
precessing warped accretion disc. The changing orientation 
of this disc warp will cause the irra diation of the secondary 
to vary (see also lVrtilek et aLlbOOSl ). 

To estimate t he effects of irradiation for Cyg X-2, we 
use the approach o f lWade fc Hornd | 19881 ). If the absorption 
centre of the secondary star is displaced from the centre of 
mass by an amount Ar = fr2, where r2 is the radius of the 
secondary star and / is the fractional displacement, then: 



AK2 



Ar 

a2 



K2 



where a2 is the separation between the secondary and the 
binary centre of mass. Based on our ne w valu e for q, a2 = 
0.70a. Using the formula of IPaczvhskil (|l97lh . the radius 
of the secondary is 0.31a. Therefore, AK2 = 0.443/7^2 • In 
the extreme case, where there is no absorption from the 
inner hemisphere of the secondary star, / = 4/(37r), and 
AK2 = 0.197^2 2i 15 km s"^ Therefore, we find that the 
measured value of K2 can be as much as ~15 km s~^ greater 
than the true value. This difference is more than enough to 
explain the different values found by various authors. 

In order to investigate the effect of irradiation on the 
absorption features of the secondary in Cyg X-2, we show 
in Fig. [TU] the EW of the Fe I A5227.15 absorption line as 
a function of the EW of He 11 A4686 (which we take as 
indicative of the intensity of the X-rays intercepted by the 
secondary star). This shows that, for this line at least, the 
stronger the He 11 A4686 emission line (and hence X-ray 
flux), the weaker the secondary absorption line. 

On the other hand, a Doppler tomogram of the 
Ca I A4227 absorption line shows that the peak occurs at 
a value consistent with that found in our RV analysis. How- 
ever, the EW of this line does not show any trend when 
plotted against the EW of the He 11 A4686 emission line. 
This may suggest that this particular line is less suscepti- 
ble to the effects of irradiation than other hues (e.g. the 
Fe I A5227.15 hne), and that the radial velocity determined 
from this line more accurately reflects the motion of the cen- 
tre of mass of the secondary star. Because irradiation of the 
secondary star causes the measured K2 to be greater than 
the true value, our measurement of a lower value may be 
closer to t he true value than th e estimates of ICowlev et al.l 
l|l979h and lCasares et all (|l998l '). 

We also note from Ta ble [4] that observations obt ained 
over a shorter time base (jCrampton fc CowlevI Il980l . and 
this work) appear to suggest systematically l ower values of 
K2 than those obtained o ver a longer period (jCowlev et al.l 
I1979I : [c asares et al.lll998l l. We speculate that observations 
obtained over a longe r timebase, and wid er range of disc 
precession phase (e.g. IClarkson et al.ll2003l ). result in more 
variable irradiation of the secondary star. It is possible that 
the changing values of AK2 that result from this may ex- 
plain the variation in K2 we observe. 



4.2 Primary mass and distance 

Our results indicate that the mass of the NS in Cyg X- 
2 {Ml — 1.5 ± 0.3 Mq) is typical of a canonical NS, a 
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Figure 10. EW of the Fe I A5227.15 absorption line vs. the EW 
of the He 11 A4686 emission line. Assuming that the EW of the 
He II A4686 emission line is correlated with the intensity of the 
X-rays intercepted by the secondary star, this relationship shows 
that the EW of the Fe I A5227.15 absorption line decreases with 
increasing X-ray flux. 



downward revision of the previou s estimates of a 95% con- 
fidence lower limi t of 1.88 M^ta (ICasares et all Il998l l and 
1.78 ± 0.23 Mq (|Orosz fc Kuulkerj 119991 )" It is interest- 
ing to note that the LMXB 2S 0921-630 has also recently 
had its primary mas s estimate reduced, from 3.2 ± 1.2 Mf^ 
llShahbaz et al.ll2004h to 1. 37±0.13 M^ llShahbaz fc WatsonI 
[^^lmdT44±0.10 M0 (|Steeghs fc Jonkeij|2007l ). Our dv- 
namical mass estimate f or the NS in Cyg X-2 is consistent 
with the mass estimate of lTitarchuk fc Shaposhnikovl (|2002l ) 
(Ml = 1.44 ± 0.06 Mq), based on observations of type-I X- 

ray bursts from this system. 

Calculations by lOrosz fc Kuulke73 (|l999l ) found that 
there are inconsistencies between the distance based on their 
optical ob servations (d — 7.2 ± 1.1 kpc) and the distance es- 
timate of ISmald i 19981 ) based on X-ray observations of a 
type-I radius-expansion burst (d = 11.6 ± 0.3 kpc). To ad- 
dress this, we now re-calculate the distance and observed 
absolute visual magnitude (Mv) of Cygnus X-2, and com- 
pare the latter with the expected value based on our new 
parameters. 

Taking t he peak X-ray flux of 1.52 x 10~* erg cm^^ s^^ 
(|Smalelll99a ) to be at the Eddington limit, and with Mi = 
1.5 ± 0.3 Mq, the Eddington luminosity of the system is 
2.0 ± 0.4 X 10''* erg s"'^ (assuming cosmic composition) and 
the distance to the source is therefore 10. 4 ± 1.0 kpc. With 
y = 1 4.8 mag, and Ay = 1.24± 0.22 mag (|McChntock et al.l 
1 19841 : [O rosz fc Kuulkerslll999l ). this gives A/v = — 1.5 ± 
0.3 mag for the system. 

Assuming a secondary star with a surface temper - 
ature of T2 = 7000 ± 250 K (|Orosz fc Kuulkerd Il999l ). 
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Ml = 1.5 ± 0.3 Mq and q = 0.42 ± 0.06, the rad ius of the 
Roch e lobe of the secondary star is 7.6 ± 0.7 Rq jEggletonl 
1 19831 ). Hence, we calculate that the Afv of the secondary 
should be —0.47 ± 0.25 mag. With a y- band disc fraction 
(fcy = /disc/(/2 + ./disc)) of 0.30 ± 0.05 jOrosz fc Kuulkerd 
1 19991 ). the absolute magnitude of the system should there- 
fore be —0.86 ± 0.26 mag. This is marginally consistent 
with the expected value for this system, calculated above 
assuming a distance of 10.4 ± 1.0 kpc. These estimates 
would better agree if the secondary temperature was slightly 
higher. Indeed, we note that the spectral classification of 
ICasares et al.l (jl998h is A9 ± 2 subtypes, and an A7 star 
(with a surface temperature of 7700 ± 300 K) would yield 
My = -0.88 ± 0.26 mag for the secondary and -1.27 ± 
0.27 mag for the system as a whole, leading to a distance of 
9.4 ± 1.5 kpc. Hence, we conclude that the absolute magni- 
tude calculated from the observed apparent magnitude and 
the distance based on X-ray measurements is consistent with 
the expected absolute magnitude based on the spectral type, 
surface temperature, disc fraction, primary mass and mass 
ratio. 

For many systems, the secondary star dominates the 
infrared light, and fitting the ellipsoidal modulation in the 
IR light curves can be used to better constrain the incli- 
nation. For Cyg X-2, such a scenario would be particu- 
larly desirable, as the modulation in the optical light curves 
is dominated by variable contamination from the accre- 
tion disc. The 2MASS A'-band magnitu de of Cyg X-2 is 
~13. 05 mag. Using E{B - V) = 0.4 jMcClintock et all 
1 19841 ). Ak/E{B -V) = 0.38 (jSavage fc Mathislll979l l aM 
a distance of 10.4 kpc, the dereddended absolute JT-band 
magnitude is AIk = — 2.2 ± 0.2 mag. If the secondary is an 
A7 star, then V — if ~ 0.8 mag. Using the absolute l/-band 
magnitude for the secondary of — 0.88±0.26 mag (calculated 
above), Mk for the secondary should be ~ — 1.7±0.3 mag. To 
achieve the observed Mk for the system of —2.2 ± 0.2 mag 
would require a 7^-band disc fraction of A:k = 0.371q'25. 
Therefore, it would appear that even in the K-hand, the ac- 
cretion disc contributes a considerable fraction of the light. 

4.3 Doppler tomography 

Our Doppler tomography shows that the peak of the He ll 
emission originates on the inner hemisphere of the donor 
star. This is consistent with the results of IVrtilek et al.l 
(|2003l ). who find that the emission lines in their UV spectra 
also appear to originate in the secondary star. There appears 
to be little or no emission from the gas stream/accr etion disc 
impact region. Curiouslv. lOrosz fc KuulkersI (|l999l ) find that 
the optical light curves show no evidence for heating of the 
secondary star, yet our He ll Doppler tomograms show that 
the secondary is certainly intercepting X-rays from the pri- 
mary. It is unclear how the secondary could be irradiated by 
the central source, and yet the light curves show no evidence 
for heating. 

The situation is reversed in the N ill A4640.64 tomo- 
gram (Fig. O , where the majority of the emission arises 
at the location of the hotspot, between the gas stream tra- 
jectory and the velocity of the accretion disc along the gas 
stream. 

Generally, the Bowen blend contains several lines, with 
the strongest normally being the N ill A4640.64 line. How- 



ever, in our data, this is the only line which is clearly visi- 
ble. The location of the emission in the Doppler tomogram 
shows that significant Bowen blend emission can occur from 
sources other than the second ary star in LMXBs. Thi s is in 
contrast to the suggestion by ICornelisse et al.l (|2008l ) that 
the velocity of the Bowen blend emission can be used uni- 
versally to determine the orbital velocity of the secondary 
star. 



5 CONCLUSIONS 

By cross-correlating our spectra against template spectra, 
we find a value for the projected radial velocity semi- 
amplitude of the se condary which is si gnificantly lower 
then that found by ICasares et al.l (|l998h . Based on this 
new value for K2, we find a larger value for the mass ra- 
tio, using the rotational broadening value calculated by 
ICasares et al.l (Il998l). Combined w ith the inclination esti- 
mates bv lOrosz fc KuulkersI (| 19991 ). this implies that the 
mass of the primary in this system is 1.5 ± 0.3 Mq, 
less massive than previously determin ed, although still 
consis tent with the mass estimate of lOrosz fc Kuulkerg 
(Il999h at the Iff level, and simil ar to the value found by 
iTitarchuk fc Shaposhnikovl (|2002l ) based on observations of 
type-I X-ray bursts. 

We have presented Doppler tomograms of Cyg X-2, the 
first Doppler tomography of such a long period system. Our 
tomography of the He 11 A4686 emission line shows that the 
majority of the emission is from near the irradiated face 
of the secondary, with much weaker emission from the ac- 
cretion disc. Doppler tomography of the N ill A4640 emis- 
sion line shows that the emission is primarily from the gas 
stream/accretion disc impact region. The N ill tomography 
demonstrates that care must be taken when interpreting 
Bowen blend tomograms, so that the hotspot is not con- 
fused with the secondary, particularly for systems where the 
ephemeris is not accurately known. 

The measured value of K2 can vary by up to 15 km s^^, 
depending on the irradiation of the secondary star by the 
central X-ray source, and this may explain the differences in 
K2 found by different authors. Further investigation of the 
effects of irradiation of the secondary star - for example, by 
performing a series of radial velocity studies over a range of 
disc precession phases - is required to fully understand how 
large this effect is on the mass estimate of the neutron star 
in Cyg X-2. 
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